Research Note: Possible influence of thermal selection on patterns of HSP70 and HSP90 gene polymorphisms in Thai indigenous and local chicken breeds and red junglefowls.

The thermal stress caused by global climate change adversely affects the welfare, productivity, and reproductive performance of farm animals, including chickens, and causes substantial economic losses. However, the understanding of the genetic basis of the indigenous chicken adaptation to high ambient temperatures is limited. Hence, to reveal the genetic basis of thermal stress adaptation in chickens, this study investigated polymorphisms in the heat shock protein 70 (HSP70) and HSP90 genes, known mechanisms of cellular defense against thermal stress in indigenous and local chicken breeds and red junglefowls in Thailand. The result revealed seven alleles of the HSP70 gene. One allele exhibited a missense mutation, where an amino acid changed from Asn to His in the substrate-binding and peptide-binding domains, which is exclusive to the Lao Pa Koi chicken breed. Twenty new alleles with silent mutations in the HSP90 gene highlighted its greater complexity. Despite this diversity, distinct population structures were not found for either HSP70 or HSP90, which suggests incomplete impact on the domestication process and selection. The low genetic diversity, shown by the sharing of alleles between red junglefowls and Thai indigenous and local chicken breeds, aligns with the hypothesis that these alleles have undergone selection in tropical regions, such as Thailand. Selection signature analysis suggests the purifying selection of HSP70 for thermotolerance. This study provides valuable insights for enhancing the conservation of genetic resources with thermotolerant traits, which are essential for developing breeding programs to increase poultry production in the context of global climate change.

Genetic insights: mapping sex-specific loci in Siamese cobra (Naja kaouthia) sheds light on the putative sex determining region.

The location of female-specific/linked loci identified in Siamese cobra (Naja kaouthia) previously has been determined through in silico chromosome mapping of the Indian cobra genome (N. naja) as a reference genome. In the present study, we used in silico chromosome mapping to identify sex-specific and linked loci in Siamese cobra. Many sex-specific and sex-linked loci were successfully mapped on the Z sex chromosome, with 227 of the 475 specific loci frequently mapped in a region covering 57 Mb and positioned at 38,992,675-95,561,177 bp of the Indian cobra genome (N. naja). This suggested the existence of a putative sex-determining region (SDR), with one specific locus (PA100000600) homologous to the TOPBP1 gene. The involvement of TOPBP1 gene may lead to abnormal synaptonemal complexes and meiotic chromosomal defects, resulting in male infertility. These findings offer valuable insights into the genetic basis and functional aspects of sex-specific traits in the Siamese cobra, which will contribute to our understanding of snake genetics and evolutionary biology.

Disclosing the hidden nucleotide sequences: a journey into DNA barcoding of raptor species in public repositories.

BACKGROUND: In nucleotide public repositories, studies discovered data errors which resulted in incorrect species identification of several accipitrid raptors considered for conservation. Mislabeling, particularly in cases of cryptic species complexes and closely related species, which were identified based on morphological characteristics, was discovered. Prioritizing accurate species labeling, morphological taxonomy, and voucher documentation is crucial to rectify spurious data. OBJECTIVE: Our study aimed to identify an effective DNA barcoding tool that accurately reflects the efficiency status of barcodes in raptor species (Accipitridae). METHODS: Barcode sequences, including 889 sequences from the mitochondrial cytochrome c oxidase I (COI) gene and 1052 sequences from cytochrome b (Cytb), from 150 raptor species within the Accipitridae family were analyzed. RESULTS: The highest percentage of intraspecific nearest neighbors from the nearest neighbor test was 88.05% for COI and 95.00% for Cytb, suggesting that the Cytb gene is a more suitable marker for accurately identifying raptor species and can serve as a standard region for DNA barcoding. In both datasets, a positive barcoding gap representing the difference between inter-and intra-specific sequence divergences was observed. For COI and Cytb, the cut-off score sequence divergences for species identification were 4.00% and 3.00%, respectively. CONCLUSION: Greater accuracy was demonstrated for the Cytb gene, making it the preferred primary DNA barcoding marker for raptors.

Characterization of the marsupial endogenous retrovirus walb with a focus on satellite DNA formation.

The walbRep megasatellite DNA found in the red-necked wallaby was formed from the walb endogenous retrovirus. Our previous PCR experiments suggested the presence of walb and absence of walbRep in the genome of the tammar wallaby, which diverged from the red-necked wallaby 2-3 Mya. The results failed to exclude the possibility that certain walbRep sequences might have remained undetected owing to variation in the primer-annealing regions; therefore, the aforementioned suggestion was not confirmed. To obtain conclusive evidence, we analyzed the structure of walb sequences drawn from the tammar wallaby genome database recently updated to a chromosome-level assembly. All walb copies existed as separate DNA segments, not constituting tandem repeats. We concluded that walbRep was formed in the red-necked wallaby lineage after its divergence from the tammar wallaby. We also confirm the presence of a walb copy with an anomalistic, complex structure and propose a simple model for its generation mechanism.

Hi-C sequencing unravels dynamic three-dimensional chromatin interactions in muntjac lineage: insights from chromosome fusions in Fea's muntjac genome.

Eukaryotes have varying numbers and structures of characteristic chromosomes across lineages or species. The evolutionary trajectory of species may have been affected by spontaneous genome rearrangements. Chromosome fusion drastically alters karyotypes. However, the mechanisms and consequences of chromosome fusions, particularly in muntjac species, are poorly understood. Recent research-based advancements in three-dimensional (3D) genomics, particularly high-throughput chromatin conformation capture (Hi-C) sequencing, have allowed for the identification of chromosome fusions and provided mechanistic insights into three muntjac species: Muntiacus muntjak, M. reevesi, and M. crinifrons. This study aimed to uncover potential genome rearrangement patterns in the threatened species Fea's muntjac (Muntiacus feae), which have not been previously examined for such characteristics. Deep Hi-C sequencing (31.42 × coverage) was performed to reveal the 3D chromatin architecture of the Fea's muntjac genome. Patterns of repeated chromosome fusions that were potentially mediated by high-abundance transposable elements were identified. Comparative Hi-C maps demonstrated linkage homology between the sex chromosomes in Fea's muntjac and autosomes in M. reevesi, indicating that fusions may have played a crucial role in the evolution of the sex chromosomes of the lineage. The species-level dynamics of topologically associated domains (TADs) suggest that TAD organization could be altered by differential chromosome interactions owing to repeated chromosome fusions. However, research on the effect of TADs on muntjac genome evolution is insufficient. This study generated Hi-C data for the Fea's muntjac, providing a genomic resource for future investigations of the evolutionary patterns of chromatin conformation at the chromosomal level.

Optimizing Microsatellite Marker Panels for Genetic Diversity and Population Genetic Studies: An Ant Colony Algorithm Approach with Polymorphic Information Content.

Microsatellites are polymorphic and cost-effective. Optimizing reduced microsatellite panels using heuristic algorithms eases budget constraints in genetic diversity and population genetic assessments. Microsatellite marker efficiency is strongly associated with its polymorphism and is quantified as the polymorphic information content (PIC). Nevertheless, marker selection cannot rely solely on PIC. In this study, the ant colony optimization (ACO) algorithm, a widely recognized optimization method, was adopted to create an enhanced selection scheme for refining microsatellite marker panels, called the PIC-ACO selection scheme. The algorithm was fine-tuned and validated using extensive datasets of chicken (Gallus gallus) and Chinese gorals (Naemorhedus griseus) from our previous studies. In contrast to basic optimization algorithms that stochastically initialize potential outputs, our selection algorithm utilizes the PIC values of markers to prime the ACO process. This increases the global solution discovery speed while reducing the likelihood of becoming trapped in local solutions. This process facilitated the acquisition of a cost-efficient and optimized microsatellite marker panel for studying genetic diversity and population genetic datasets. The established microsatellite efficiency metrics such as PIC, allele richness, and heterozygosity were correlated with the actual effectiveness of the microsatellite marker panel. This approach could substantially reduce budgetary barriers to population genetic assessments, breeding, and conservation programs.

Overcoming taxonomic challenges in DNA barcoding for improvement of identification and preservation of clariid catfish species.

DNA barcoding without assessing reliability and validity causes taxonomic errors of species identification, which is responsible for disruptions of their conservation and aquaculture industry. Although DNA barcoding facilitates molecular identification and phylogenetic analysis of species, its availability in clariid catfish lineage remains uncertain. In this study, DNA barcoding was developed and validated for clariid catfish. 2,970 barcode sequences from mitochondrial cytochrome c oxidase I (COI) and cytochrome b (Cytb) genes and D-loop sequences were analyzed for 37 clariid catfish species. The highest intraspecific nearest neighbor distances were 85.47%, 98.03%, and 89.10% for COI, Cytb, and D-loop sequences, respectively. This suggests that the Cytb gene is the most appropriate for identifying clariid catfish and can serve as a standard region for DNA barcoding. A positive barcoding gap between interspecific and intraspecific sequence divergence was observed in the Cytb dataset but not in the COI and D-loop datasets. Intraspecific variation was typically less than 4.4%, whereas interspecific variation was generally more than 66.9%. However, a species complex was detected in walking catfish and significant intraspecific sequence divergence was observed in North African catfish. These findings suggest the need to focus on developing a DNA barcoding system for classifying clariid catfish properly and to validate its efficacy for a wider range of clariid catfish. With an enriched database of multiple sequences from a target species and its genus, species identification can be more accurate and biodiversity assessment of the species can be facilitated.

MicrosatNavigator: exploring nonrandom distribution and lineage-specificity of microsatellite repeat motifs on vertebrate sex chromosomes across 186 whole genomes.

Microsatellites are short tandem DNA repeats, ubiquitous in genomes. They are believed to be under selection pressure, considering their high distribution and abundance beyond chance or random accumulation. However, limited analysis of microsatellites in single taxonomic groups makes it challenging to understand their evolutionary significance across taxonomic boundaries. Despite abundant genomic information, microsatellites have been studied in limited contexts and within a few species, warranting an unbiased examination of their genome-wide distribution in distinct versus closely related-clades. Large-scale comparisons have revealed relevant trends, especially in vertebrates. Here, "MicrosatNavigator", a new tool that allows quick and reliable investigation of perfect microsatellites in DNA sequences, was developed. This tool can identify microsatellites across the entire genome sequences. Using this tool, microsatellite repeat motifs were identified in the genome sequences of 186 vertebrates. A significant positive correlation was noted between the abundance, density, length, and GC bias of microsatellites and specific lineages. The (AC)n motif is the most prevalent in vertebrate genomes, showing distinct patterns in closely related species. Longer microsatellites were observed on sex chromosomes in birds and mammals but not on autosomes. Microsatellites on sex chromosomes of non-fish vertebrates have the lowest GC content, whereas high-GC microsatellites (≥ 50 M% GC) are preferred in bony and cartilaginous fishes. Thus, similar selective forces and mutational processes may constrain GC-rich microsatellites to different clades. These findings should facilitate investigations into the roles of microsatellites in sex chromosome differentiation and provide candidate microsatellites for functional analysis across the vertebrate evolutionary spectrum.

Kangaroo endogenous retrovirus (KERV) forms megasatellite DNA with a simple repetition pattern in which the provirus structure is retained.

The kangaroo endogenous retrovirus (KERV) was previously reported to have undergone a rapid copy number increase in the red-necked wallaby; however, the mode of amplification was left to be clarified. The present study revealed that the long terminal repeat (LTR) (0.6 kb) and internal region (2.0 kb) of a provirus are repeated alternately, forming megasatellite DNA which we named kervRep. This repetition pattern was the same as that observed for walbRep, megasatellite DNA originating from another endogenous retrovirus. Their formation process can be explained using a simple model: pairing slippage followed by homologous recombination. This model features that the initial step is triggered by the presence of two identical sequences within a short distance; the possession of LTRs by endogenous retroviruses fulfills this condition. The discovery of two cases suggests that formation of this type of satellite DNA is one of non-negligible effects of endogenous retroviruses on their host genomes.

Marsupial genome analysis suggests that satellite DNA formation from walb endogenous retrovirus is an event specific to the red-necked wallaby.

We recently identified walbRep, a satellite DNA residing in the genome of the red-necked wallaby Notamacropus rufogriseus. It originates from the walb endogenous retrovirus and is organized in a manner in which the provirus structure is retained. The walbRep repeat units feature an average pairwise nucleotide identity as high as 99.5%, raising the possibility of a recent origin. The tammar wallaby N. eugenii is a species estimated to have diverged from the red-necked wallaby 2-3 million years ago. In PCR analyses of these two and other related species, walbRep-specific fragment amplification was observed only in the red-necked wallaby. Sequence database searches for the tammar wallaby resulted in sequence alignment lists that were sufficiently powerful to exclude the possibility of walbRep existence. These results suggested that the walbRep formation occurred in the red-necked wallaby lineage after its divergence from the tammar wallaby lineage, thus in a time span of maximum 3 million years.

Marsupial satellite DNA as faithful reflections of long-terminal repeat retroelement structure.

Long terminal repeat (LTR) retroelements, including endogenous retroviruses, are one of the origins of satellite DNAs. However, the vast majority of satellite DNAs originating from LTR retroelements consists of parts of the element. In addition, they frequently contain sequences unrelated to that element. Here we report a novel marsupial satellite DNA (named walbRep) that contains, and consists solely of, the entire sequence of an LTR retroelement (the walb element). As is common with LTR retroelements, walb copies exhibit length variation. We focused on the abundance of copies of a specific length (2.7 kb) in the genome of the red-necked wallaby. Cloning and analyses of long genomic DNA fragments revealed a satellite DNA in which the LTR sequence (0.4 kb) and the sequence of the internal region of a nonautonomous walb copy (2.3 kb) were repeated alternately. The junctions between these two components exhibited the same end-to-end arrangements as those in the walb element. This satellite organization could be accounted for by a simple formation model that includes slippage during chromosome pairing followed by homologous recombination but does not invoke any other types of rearrangements. We discuss the possible reasons why satellite DNAs having such structures are rarely found in mammals.

An endogenous retrovirus presumed to have been endogenized or relocated recently in a marsupial, the red-necked wallaby.

An albino infant wallaby was born to a mother with wild-type body color. PCR and sequencing analyses of TYR (encoding tyrosinase, which is essential for melanin biosynthesis) of this albino wallaby revealed a 7.1-kb-long DNA fragment inserted in the first exon. Since the fragment carried long terminal repeats, we assumed it to be a copy of an endogenous retrovirus, which we named walb. We cloned other walb copies residing in the genomes of this species and of another wallaby species. The copies exhibited length variation, and the longest copy (>8.0 kb) contained open reading frames whose deduced amino acid sequences were well aligned with those of gag, pol, and env of retroviruses. It is unknown through which of the following likely processes the walb copy was inserted into TYR: endogenization (infection of a germline cell by an exogenous virus), reinfection (infection by a virus produced from a previously endogenized provirus), or retrotransposition (intracellular relocation of a provirus). In any case, the insertion into TYR is considered to have been a recent event on an evolutionary timescale because albino mutant alleles generally do not persist for long because of their deleterious effects in wild circumstances.

The medaka fish Tol2 transposable element is in an early stage of decay: identification of a nonautonomous copy.

The majority of DNA-based transposable elements comprise autonomous and nonautonomous copies, or only nonautonomous copies, where the autonomous copy contains an intact gene for a transposase protein and the nonautonomous copy does not. Even if autonomous copies coexist, they are generally less frequent. The Tol2 element of medaka fish is one of the few elements for which a nonautonomous copy has not yet been found. Here, we report the presence of a nonautonomous Tol2 copy that was identified by surveying the medaka genome sequence database. This copy contained three local sequence alterations that affected the deduced amino acid sequence of the transposase: a deletion of 15 nucleotides resulting in a deletion of 5 amino acids, a base substitution causing a single amino acid change, and another base substitution giving rise to a stop codon. Transposition assays using cultured human cells revealed that transposase activity was reduced by the 15-nucleotide deletion and abolished by the nonsense mutation. This is the first example of a nonautonomous Tol2 copy. Thus, Tol2 is in an early stage of decay in the medaka genome, and is therefore a unique element to observe an almost complete decay process that progresses in natural populations.

Replacement of owl monkey centromere satellite by a newly evolved variant was a recent and rapid process.

Alpha satellite DNA is a major DNA component of primate centromeres. We previously reported that Azara's owl monkey has two types of alpha satellite DNA, OwlAlp1 and OwlAlp2. OwlAlp2 (344 bp) exhibits a sequence similarity throughout its entire length with alpha satellite DNA of closely related species. OwlAlp1 (185 bp) corresponds to the part of OwlAlp2. Based on the observation that the CENP-A protein binds to OwlAlp1, we proposed that OwlAlp1 is a relatively new repetitive DNA that replaced OwlAlp2 as the centromeric satellite DNA. However, a detailed picture of the evolutionary process of this centromere DNA replacement remains largely unknown. Here, we performed a phylogenetic analysis of OwlAlp1 and OwlAlp2 sequences, and also compared our results to alpha satellite DNA sequences of other primate species. We found that: (i) OwlAlp1 exhibits a higher similarity to OwlAlp2 than to alpha satellite DNA of other species, (ii) OwlAlp1 has a single origin, and (iii) sequence variation is lower in OwlAlp1 than in OwlAlp2. We conclude that OwlAlp1 underwent a recent and rapid expansion in the owl monkey lineage. This centromere DNA replacement could have been facilitated by the heterochromatin reorganization that is associated with the adaptation of owl monkeys to a nocturnal lifestyle.

A mutant gene for albino body color is widespread in natural populations of tanuki (Japanese raccoon dog).

Albino mutants (white coat and red eyes) of tanuki (Nyctereutes procyonoides viverrinus) have been repeatedly found in the Central Alps area of Japan. We recently reported that an albino tanuki from Iida, a city in this area, lacks the third exon of the TYR gene encoding tyrosinase, which is essential for melanin synthesis. The absence of this exon was due to the chromosomal deletion of a complex structure. In the present study, we analyzed TYR of another albino tanuki that was found in Matsusaka, a city located outside the mountainous area. In this animal, the third exon was also lost, and the loss was due to a deletion in which the structure was identical to that of the Iida mutant. Our results indicate, in consideration of the complex structure of the deletion, that the two albino animals inherited a single deletion that arose in their common ancestor. Iida and Matsusaka are approximately 170 km apart. This is, to our knowledge, the first report of an albino mutant gene that is widely distributed in mammalian natural populations. As the origin of this mutation is not known, the distance covered by the mutant gene remains unclear. If we assume that the mutation occurred halfway between Iida and Matsusaka, we can predict the migration distance to be approximately 85 km; however, if the mutation occurred at any other place, a longer distance would be predicted. Natural selection against albino tanuki may be relaxed because of a recent increase in food resources and refuge in urban areas.

The Heterochromatin Block That Functions as a Rod Cell Microlens in Owl Monkeys Formed within a 15-Myr Time Span.

In rod cells of many nocturnal mammals, heterochromatin localizes to the central region of the nucleus and serves as a lens to send light efficiently to the photoreceptor region. The genus Aotus (owl monkeys) is commonly considered to have undergone a shift from diurnal to nocturnal lifestyle. We recently demonstrated that rod cells of the Aotus species Aotus azarae possess a heterochromatin block at the center of its nucleus. The purpose of the present study was to estimate the time span in which the formation of the heterochromatin block took place. We performed three-dimensional hybridization analysis of the rod cell of another species, Aotus lemurinus. This analysis revealed the presence of a heterochromatin block that consisted of the same DNA components as those in A. azarae. These results indicate that the formation was complete at or before the separation of the two species. Based on the commonly accepted evolutionary history of New World monkeys and specifically of owl monkeys, the time span for the entire formation process was estimated to be 15 Myr at most.

Complex intragene deletion leads to oculocutaneous albinism in tanuki (Japanese raccoon dog).

Tanuki (Nyctereutes procyonoides viverrinus), or Japanese raccoon dog, is a canine native to Japan. Tanuki with complete oculocutaneous albinism are relatively frequent in mountainous areas of mainland Japan. Tyrosinase, which is encoded by the TYR gene, is an enzyme essential for the biosynthesis of melanin pigment. We examined the structure and nucleotide sequence of TYR in an albino tanuki and found that the third exon was removed due to a deletion of approximately 11 kb. In addition, two nonsynonymous nucleotide substitutions were found in the fifth exon. These mutations are possible causes of the albino phenotype; however, the order of occurrence is unclear. Even if the 11-kb deletion was not the first of these mutations, it is considered to cause a total loss of the tyrosinase function because the third exon carries codons for one of the two copper-binding sites of tyrosinase and these sites are essential for the enzyme function. Intriguingly, the deletion was not a simple removal of an 11-kb segment: an internal portion was retained as a segment in the reverse orientation. We propose possible formation processes for this mutation that involve multiple DNA scission events, or an inversion followed by a deletion.

Baboon bearing resemblance in pigmentation pattern to Siamese cat carries a missense mutation in the tyrosinase gene.

An infant hamadryas baboon exhibiting an albino phenotype-white body hair and red eyes-was born to parents with wild-type body color. Pigmentation on some parts of its body surfaced during childhood and progressed with age. This baboon in adulthood has gray hair on parts of its body, such as the tail, distal portion of the legs, and face, with the remainder being white. This pigmentation pattern resembles that of the Siamese cat and the Himalayan variants of the mouse and the mink. The distinguishing phenotypes in these animals are known to be caused by a temperature-sensitive activity of tyrosinase, an enzyme essential for biosynthesis of melanin. We sequenced all the five exons of the tyrosinase (TYR) gene of this albino baboon, which were amplified by PCR, and found a base substitution leading to alteration of the 365th amino acid from Ala to Thr. Tyrosinase requires copper as a cofactor for its enzyme function. It has two copper-binding sites, the second of which contains His residues in positions 363 and 367 that are critical to its function. Thus, p.(Ala365Thr) due to a mutation in the TYR gene is a likely candidate for the cause of the albino phenotype in this baboon.

The CENP-B box, a nucleotide motif involved in centromere formation, has multiple origins in New World monkeys.

Centromere protein B (CENP-B), a protein participating in centromere formation, binds to centromere satellite DNA by recognizing a 17-bp motif called the CENP-B box. This motif is found in hominids (humans and great apes) at an identical location in repeat units of their centromere satellite DNA. We have recently reported that the CENP-B box exists at diverse locations in three New World monkey species (marmoset, squirrel monkey and tamarin). However, the evolutionary origin of the CENP-B box in these species was not determined. It could have been present in a common ancestor, or emerged multiple times in different lineages. Here we present results of a phylogenetic analysis of centromere satellite DNA that support the multiple emergence hypothesis. Repeat units almost invariably formed monophyletic groups in each species and the CENP-B box location was unique for each species. The CENP-B box is not essential for the immediate survival of its host organism. On the other hand, it is known to be required for de novo centromere assembly. Our results suggest that the CENP-B box confers a long-term selective advantage. For example, it may play a pivotal role when a centromere is accidentally lost or impaired.

Structural variations of subterminal satellite blocks and their source mechanisms as inferred from the meiotic configurations of chimpanzee chromosome termini.

African great apes have large constitutive heterochromatin (C-band) blocks in subtelomeric regions of the majority of their chromosomes, but humans lack these. Additionally, the chimpanzee meiotic cell division process demonstrates unique partial terminal associations in the first meiotic prophase (pachytene). These are likely formed as a result of interaction among subtelomeric C-band blocks. We thus conducted an extensive study to define the features in the subtelomeric heterochromatic regions of chimpanzee chromosomes undergoing mitotic metaphase and meiotic cell division. Molecular cytogenetic analyses with probes of both subterminal satellite DNA (a main component of C-band) and rDNA demonstrated principles of interaction among DNA arrays. The results suggest that homologous and ectopic recombination through persistent subtelomeric associations (post-bouquet association observed in 32% of spermatocytes in the pachytene stage) appears to create variability in heterochromatin patterns and simultaneously restrain subtelomeric genome polymorphisms. That is, the meeting of non-homologous chromosome termini sets the stage for ectopic pairing which, in turn, is the mechanism for generating variability and genomic dispersion of subtelomeric C-band blocks through a system of concerted evolution. Comparison between the present study and previous reports indicated that the chromosomal distribution rate of sutelomeric regions seems to have antagonistic correlation with arm numbers holding subterminal satellite blocks in humans, chimpanzees, and gorillas. That is, the increase of subterminal satellite blocks probably reduces genomic diversity in the subtelomeric regions. The acquisition vs. loss of the subtelomeric C-band blocks is postulated as the underlying engine of this chromosomal differentiation yielded by meiotic chromosomal interaction.